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We construct the boundary conformal field theory that describes the low-temperature behavior of 
the two-channel Anderson impurity model. The presence of an exactly marginal operator is shown 
to generate a line of stable fixed points parameterized by the charge valence n c of the impurity. We 
calculate the exact zero-temperature entropy and impurity thermodynamics along the fixed line. 
We also derive the critical exponents of the characteristic Fermi edge singularities caused by time- 
dependent hybridization between conduction electrons and impurity. Our results suggest that in 
the mixed-valent regime (n c 7^ 0, 1) the electrons participate in two competing processes, leading to 
frustrated screening of spin and channel degrees of freedom. By combining the boundary conformal 
field theory with the Bethe Ansatz solution we obtain a complete description of the low-energy 
dynamics of the model. 
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I. INTRODUCTION 

In recent years, a growing class of materials has been 
shown to exhibit metallic behaviors that violate Lan- 
dau's Fermi liquid theory. 1 Examples include the "nor- 
mal phases" of underdoped and optimally doped high- 
T c cuprates^ a variety of (quasi) one-dimensional con- 
ductors - ranging from single-wall carbon nanotubesSi 
to the Bechgaard salts^ - several artificially designed 
nanostructuresp as well as certain cerium- and uranium- 
based heavy fermion alloys > e In light of the spectacular 
success of Landau's theory in explaining the properties 
of conventional metals, the proliferation of experimental 
systems that depart from its predictions presents a chal- 
lenge to the theorist. Not forming a Fermi liquid, the 
mobile electrons of these systems cannot be adiabatically 
connected to a noninteracting electron gas, and their de- 
scription requires a theoretical approach of a different 
brand. 

A particularly intriguing case is that of the heavy 
fermion compound UBei3. Like several other metallic 
materials containing rare earth or actinide ions with par- 
tially filled f-electron shells UBei3 exhibits manifest non- 
Fermi liquid behavior at low temperatures. In particular, 
the specific heat shows a T In T behavior all the way down 
to the superconducting transition (at roughly 1K)jI sug- 
gestive of a two-channel overscreened Kondo effect driven 
by the f-electrons of the uranium ions. There are two 
possible integer valence configurations for uranium em- 
bedded in a Bei3 host: a Tq magnetic spin doublet in 
the 5/ 3 configuration and a F3 electrical quadrupolar 
doublet in the 5/ 2 configuration, the two levels being 
separated by an energy e. Arguing that the extremely 
weak magnetic field dependence of the material^ excludes 



the usual magnetic Kondo effect, Goy& proposed that 
the observed anomalies instead derive from a quench- 
ing of the quadrupolar degrees of freedom by the lo- 
cal orbital motion of the Y% conduction electrons in the 
material (with the electron spin a =],{ providing the 
two "channels" required for overscreening) . The pro- 
posal has been controversial,?' however, and experimen- 
tal data on the nonlinear magnetic susceptibility indicate 
that the low-lying magnetic excitations are rather pre- 
dominantly dipolar in character— Whether the energy 
splitting e between the two doublets is sufficiently large 
for a quadrupolar Kondo scenario to become viable is an 
additional unresolved issue. Aliev et a/— have suggested 
that the magnetic and quadrupolar doublets may in fact 
be near degenerate, leading instead to a mixed-valent 
state with a novel type of interplay between quadrupo- 
lar and magnetic Kondo-type screening. Indirect support 
for this hypothesis comes from non- linear susceptibility 11 
and other— measurements on the thoriated compound 
Uo.gTho.iBei3. The presence of thorium ions suppresses 
the superconducting transition in the undoped material 
and one observes instead a crossover from a magnetic to 
a quadrupolar groundstate at low temperatures. 

The simplest model that captures a possible mixed- 
valent regime with both magnetic and quadrupolar 
character is the two-channel Anderson single-impurity 
model— iiii^i In this model the conduction electrons 
carry both spin (o =] , [) and quadrupolar (a = ±) quan- 
tum numbers, and hybridize via a matrix element V with 
a local uranium ion. The ion is modeled by a quadrupo- 
lar (5/ 2 ) doublet created by a boson operator b' a and a 
magnetic (5/ 3 ) doublet created by a fermion operator /J. 
Strong Coulomb repulsion implies that the localized /- 
levels can carry at most one electron, and this condition 
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is implemented as an operator identity for the pseudopar- 
ticles: f\$ a + b' ! a b a — 1, where a denotes the conjugate 
representation. With these provisos the Hamiltonian is 
written as 

H = Hbulk + H ion + Hhybr , (1) 

where 

Hbulk = J dx :^l a (x)(-id x )ip aa (x): (2) 
H ion = e s f^f a + e q b a b s (3) 

H hybr = ^[V>L(0)4/a + /<M«r(0)] ■ (4) 

Here the conduction electrons are described by one- 
dimensional fields il)< al7 (x), matching the T% representa- 
tion (appearing in the reduction of T3 x Tq and thus 
coupling to the impurity). The electron fields are chiral, 
obtained - as usuali& - by reflecting the outgoing radial 
(x > 0), r 8 -component of the original three-dimensional 
problem to the negative x-axis and disregarding com- 
ponents that do not couple to the ion at x = 0. The 
spectrum is linearized around the Fermi level and the 
Fermi velocity is set to unity with the resulting density 
of states being p = 1/ (2tt) . The normal ordering is taken 
with respect to the filled Fermi sea, and the energies e s 
and e q are those of the magnetic and quadrupolar dou- 
blets respectively. 

Recently two of us presented an exact solution of the 
model, based on a Bethe Ansatz construction »il A com- 
plete determination of the energy spectrum and the ther- 
modynamics was given, allowing a full description of the 
evolution of the impurity from its high temperature be- 
havior with all four impurity states being equally pop- 
ulated down to the low energy dynamics characterized 
by a line of fixed point Hamiltonians H*(s,T) where 
e = e s — s q and T — npV 2 (T is held fixed in what 
follows) . 

It is found that the line of fixed points is character- 
ized by a zero-temperature entropy Sj^p = fcs lnv2 and 
a specific heat C™ p ~ TlnT typical of the two-channel 
Kondo fixed point. However, the physics along the line 
varies with e. Consider n c — (f^f a ), the amount of 
charge localized at the impurity. For e < p, — T, one 
finds that n c « 1, signaling the magnetic integral va- 
lence regime. At intermediate temperatures a magnetic 
moment forms which undergoes frustrated screening as 
the temperature is lowered, leading to zero-temperature 
anomalous entropy and anomalous specific heat. For 
e > fi+T the system is in the quadrupolar integral regime 
and a quadrupolar moment forms. In the mixed valence 
regime, |e — ji\ < T, similar low-temperature behavior is 
observed though without the intermediate regime of mo- 
ment formation. In more detail, each point on the line of 
fixed points is characterized by two energy scales Ti t h{s). 
These scales describe the quenching of the entropy as the 
temperature is lowered: the first stage taking place at the 
high- temperature scale T^, quenching the entropy from 



fee In 4 to fee In 2, the second stage at 7], quenching it 
from ks In 2 to fee In a/2- In the integral valence regime, 
|e — fi\ ^> r, the two scales are well separated and as long 
as the temperature falls between these values a moment 
is present (magnetic or quadrupolar depending on the 
sign of e — fi) , manifested by a finite temperature plateau 
<Si m p = fcsln2 in the entropy. It is quenched when the 
temperature is lowered below TJ, with Ti — > Tk in this 
regime. For e = p the scales are equal, Ti(fi) = Th(ji), 
and the quenching occurs in a single stage. 

In the present paper we shall concentrate on the low- 
energy regime and give a detailed description of the line of 
fixed points in terms of boundary conformal field theory 

— the low-energy effective Hamiltonian. It was argued 
by Affleck and Ludwig 18,1 ? that the low-energy regime of 
a quantum impurity problem is described by a boundary 
conformal field theory (BCFT) with a conformally in- 
variant boundary condition replacing the dynamical im- 
purity. However, to identify the effective low-energy the- 
ory characterizing a given microscopic impurity model, 
such as Hamiltonian Q, it is necessary to employ more 
powerful methods, valid over the full energy range, that 
can connect the microscopic Hamiltonian to the effective 
low-energy theory. One needs in principle to carry out a 
full renormalization group calculation or, when available, 
use an exact solution to extract the BCFT. 

Although applicable only in the neighborhood of a low- 
temperature fixed point, a BCFT formulation has certain 
advantages. First, it provides an elegant scheme in which 
to analyze the approach to criticality. The leading scal- 
ing operators can be identified explicitly, yielding access 
to the low-temperature thermodynamics in closed ana- 
lytical form. These can be matched with the expressions 
from the Bethe- Ansatz solution allowing the determina- 
tion of the parameters appearing in the BCFT. This will 
be carried out below. Secondly — and more importantly 

— having identified the scaling operators, the asymptotic 
dynamical properties (Green's functions, resistivities, op- 
tical conductivities, etc.) can in principle be calculated. 
These results will be presented in a subsequent work. 

The paper is organized as follows. In Sec. II we re- 
view the basics of BCFT, with particular focus on ap- 
plications to quantum impurity problems. In Sec. Ill 
we derive the specific BCFT which describes the low- 
temperature physics of the two-channel Anderson model 
in Eq. Some important features of this BCFT are 

discussed, and in this section we also apply it to the Fermi 
edge singularity problem for this model. In Sec. IV we 
then combine the results obtained with the Bethe Ansatz 
solution 17 to analytically calculate the zero-temperature 
impurity entropy, the impurity contributions to the low- 
temperature specific heat, and the linear and non-linear 
impurity magnetic susceptibilities. Section V, finally, 
contains a summary and a discussion of our results. 
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II. BOUNDARY CONFORMAL FIELD THEORY 
APPROACH 

The BCFT approach to quantum impurity problems is 
well reviewed in the Iiteraturef22i2ii22i and we here only 
collect some basic results so as to fix conventions and 
notation. 

The key idea is to trade a local impurity-electron in- 
teraction for a scale invariant boundary condition on the 
critical theory that represents the extended (electron) de- 
grees of freedom. In the presence of a boundary the left- 
and right-moving parts of the fields f(t, x) in the criti- 
cal theory are identified via analytic continuation beyond 
the boundary x = 0, such that 

ip(t,x) = ip L (t,x)ip R (t,x) -> ip L (t,x)(p L (t,-x). (5) 

Thus, the boundary effectively turns the fields nonlocal. 
The consequence of this can be delineated via the oper- 
ator product expansion 

H>l{*)<Pl{?) - E {z _ z %a-a 3 ( ? (0, r) , 

z -> z* . (6) 

We have here introduced a complex variable z = t + ix, 
with t the Euclidean time and x the spatial coordinate. 
The function gj (= or 1) selects the boundary opera- 
(i) 

tors 0£ associated with the given boundary condition, 

Aj (A) being the dimension of Ojf) (fL)- It follows from 
Eqs. (0 and that an n-point function of a field <p close 
to the boundary turns into a linear combination of n- 
point functions of the chiral boundary operators . In 
particular, this implies that the boundary scaling dimen- 
sion Abound of <f that governs its autocorrelation function 
close to the boundary, 

(f(t,x)f(0,x)) - (f>(t,x))(f(0,x)) ~ i- 2 ^—" 

\t\»\x\ (7) 

is precisely given by the scaling dimension of the leading 
boundary operator appearing in Eq. ijSJ. 

This sets the strategy for treating a quantum impurity 
problem: identify first the particular boundary condition 
that plays the role of the impurity interaction. Given 
that this boundary condition is indeed scale invariant 
(together with the original bulk theory), one can then 
use the machinery of BCFT to extract the corresponding 
boundary scaling dimensions. As these determine the 
asymptotic autocorrelation functions [cf. Eq. Q] the 
finite-temperature properties due to the presence of the 
boundary (alias the impurity) are easily accessed from 
standard finite-size scaling by treating the (Euclidean) 
time as an inverse temperature. 

To identify the "right" boundary condition it is con- 
venient to exploit a well-known resulli^ from conformal 



field theory, relating the energy levels in a finite ge- 
ometry to the (boundary) scaling dimensions of opera- 
tors in the (semi-) infinite plane. More precisely, con- 
sider a conformally invariant theory defined on the strip 
{w = u + iv | — oo < u < oo, $J v ^ £} with u the 
Euclidean time and v the space coordinate. Then impose 
a conformally invariant boundary condition, call it A, at 
the edges v = and v = t, and map the strip onto the 
semi- infinite plane {z — r + ix\x ^ 0} using the confor- 
mal transformation z = exjp(irw/l) (implying boundary 
condition A at x = 0). With E the ground-state energy, 
one has 

E n = Eq H j?-, (8) 

where {E n } is the spectrum of excited energy levels in 
^ v I and {A„} is the spectrum of boundary scaling 
dimensions in the semi-infinite plane. 

The problem is thus reduced to determining the finite- 
size spectrum of the theory. In certain privileged cases 
this can be done by direct calculation. Alternatively, 
one hypothesizes a boundary condition and compares 
the consequences with known answers. An example is 
the fusion hypothesis of Affleck and Ludwig^S One here 
starts with some known, trivial boundary condition on 
the critical theory, with no coupling to the impurity. For 
free electrons carrying charge, spin (and maybe some ad- 
ditional "flavor" degrees of freedom) the spectrum orga- 
nizes into a conformal tower structure^ with the charge, 
spin and flavor sectors "glued" together so as to correctly 
represent the electrons (Fermi liquid gluing condition). 
Then, turning on the electron-impurity interaction, its 
only effect — according to the fusion hypothesis — is to 
replace this gluing condition by some new nontrivial glu- 
ing of the conformal towers. In the case of the m-channel 
Kondo problem, the new gluing condition is obtained via 
Kac-Moody fusion 25 with the spin-1/2 primary operator 
that "corresponds" to the Kondo impurity. As a con- 
sequence, the conformal tower with spin quantum num- 
ber j s is mapped onto new towers labeled by j' s , where 
j' s = \js - 1 1 , 1 3s - || + 1, • •• , min{ j s + i , m - j s - \ } . This 
is the essence of the BCFT approach, as applied to the 
Kondo problem. 

For the two-channel Anderson model in Eq. it is 
a priori less obvious which operator to use for fusion. 
Here we shall instead identify the finite-size spectrum by 
studying certain known limiting cases, guided by the ex- 
act Bethe Ansatz solution of the models As it turns 
out, the spectrum thus obtained can be reconstructed by 
fusion with the leading flavor boundary operator (repre- 
senting the channel, or quadrupolar, degrees of freedom), 
but in addition there occurs an effective renormalization 
of the charge sector. This signals the novel aspect of the 
present problem. A technical remark may here be appro- 
priate: Since the boundary scaling dimensions in Eq. (JSJ) 
are connected to energy levels in a strip with the same 
boundary condition at the two edges, we are effectively 
considering a finite-size energy spectrum with two quan- 
tum impurities present, one at each edge of the strip. 
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Formally, this can be taken care of by performing Kac- 
Moody fusion twice with the relevant primary operator 
(double fusion). This is an important point, not always 
fully appreciated in the literature. 

Let us finally mention that there exists another, more 
fundamental description of BCFTj2& based on the no- 
tion of boundary states, particularly useful when studying 
zero-temperature properties of a quantum impurity prob- 
lem. We shall give a brief exposition of it in Sec. (IV. A), 
where we use it for analyzing the zero-temperature en- 
tropy contributed by the impurity. 



III. FINITE-SIZE SPECTRUM AND SCALING 
OPERATORS 

As we discussed in the previous section, the scaling 
dimensions of the boundary operators that govern the 
critical behavior are in one-to-one correspondence with 
the levels of the finite-size spectrum of the model. In the 
present case the spectrum can in principle be constructed 
from the exact Bethe Ansatz solution in Ref. 0, but 
this requires an elaborate analysis. The reason is that 
the solution takes the form of a so called "string solu- 
tion" even in the groundstate. Since a string solution is 
valid only in the thermodynamic limit, an application of 
standard finite-size techniques^ would lead to incorrect 
results. One may instead pursue another strategy and 
identify the finite size spectrum (hence the dimensions of 
the scaling operators) that reproduces the results of the 
exact solution in the low-energy limit. In our case the 
search for the wanted spectrum can be carried out effi- 
ciently by combining symmetry arguments with known 
results for the finite-size spectra for two limiting cases of 
the Hamiltonian qTp : the ordinary Anderson model and 
the two-channel Kondo model. We shall see that the low 
energy spectrum thus obtained, Eq. (|22[1 below, leads in- 
deed to a line of fixed points characterized by entropy 
S = ks In \/2 and susceptibility Cy ~ TlnT as required 
by the Bethe Ansatz solution. The fitting of the proposed 
spectrum to the solution also allows the determination of 
all thcrmodynamically relevant parameters in the BCFT. 



A. Identifying the Critical Theory 

Let us begin by reviewing the single channel Ander- 
son and Kondo models. The latter is obtained in the 



integer- valence limit n c — > 1 (with n c = (f^f a ) mea- 
suring the charge localized at the impurity site) where 
a magnetic moment forms, signaling the entrance to the 
Kondo regime^ It is instructive to consider this case in 
some detail before proceeding to an analysis of the full 
model in Eq. 

As is well known, the ordinary Anderson models 
with integer valence n c = 1 can be mapped onto 
the single-channel Kondo model via a Schrieffer- Wolff 
transformation. 28 To 0(\/t), with £ the length of the sys- 
tem, the finite-size energy spectrum of the Kondo model 
takes the formic 



E — Eq 



+ N C ) + 



js(js + 1) 



(9) 

where the Fermi velocity has been set to unity, as in Eq. 
||TJ. Here Q € Z and j s = 0,1/2 are charge and spin 
quantum numbers - defined with respect to the filled 
Fermi sea - and labeling the corresponding Kac-Moody 
primary states*^ The positive integers N c and N s index 
the descendant levels in the associated U(l) and SU(2)i 
conformal towers. The quantum numbers are constrained 
to appear in the combinations 



Q = mod 2, 
Q = 1 mod 2, 



(10) 



which define the Kondo gluing conditions for charge and 
spin conformal towers. 

We can rewrite this result in a form immediately gener- 
alizable to the Anderson model, viewing it as the n c — > 1 
limit of the Anderson model spectrum. Redefine the 
charge quantum number in Eqs. JHJl and (|10fl : Q — > Q— L 
The offset mirrors the fact that the Anderson impurity 
carries charge (in contrast to a Kondo impurity) and, 
therefore, shifts the net amount of charge in the ground- 
state compared to the Kondo case. We thus write for the 
finite-size spectrum of the Anderson model in the n c — > 1 
limit: 



E — Eq 



:(Q-iy-- + N c + 



Jstis + 1) 



(11) 



with 



Note that the constraints in Eqs. I|12|) are the same as 



Q = mod 2, j s = 



= 1 mod 2, ] s 



1 



(12) 



5 



those for the charge and spin conformal towers of free 
electrons (Fermi liquid gluing conditions). Also note that 
we have renormalized the groundstate energy, so that 
Eq = E(Q = 0), by subtracting off a constant tt/M. 
The redefinition of charge quantum numbers implies a 
relabeling of levels: the Q-level in Eq. (constrained 
by Eqs. ljTU|) 1 corresponds to the (Q— l)-level in Eq. (fTT|) 
(with the constraint tyl'2\i ). modulo the ir/4£ shift. 

Let us now consider the spectrum of the Anderson 
model away from the Kondo limit, n c ^ 1. Assuming 

I 



E = E + - 



with Q and j s constrained by Eqs. (|12(l . Note that this 
result has been extracted from the Kondo finite-size spec- 
trum in Eqs. © and (|ll)fl simply by keeping track on the 
offset of charge with respect to the filled Fermi sea as the 
impurity acquires a nonzero charge valence n c . The addi- 
tional assumption that no relevant operator is generated 
as we move away from the Kondo limit n c = 1 (and that 
therefore the gluing conditions in Eqs. (|10f> remain un- 

I 



E = E + - 



with Fermi liquid gluing (|12(l of Q and j s . The phase shift 
5f = So-(cf) appearing in (|14H is that of an electron with 
spin a at the Fermi level, related to the average charge 
n c = at the impurity site by the Friedel-Langreth 

sum rule^ 

^F=2 n c- (15) 

Inserting Eq. I|15[) into Eq. l|14f> . we immediately recover 
our result l|13|) . 

Before returning to the full two-channel Anderson 
model, let us point out that the duality between the 
charge shift Q <-> Q — 1 in the integer valence limit n c = 1 
and the change of gluing conditions i|12fl «-> (I10|) simply 
reflects the well-known fact that the ordinary Kondo ef- 
fect can formally be described as a local potential scat- 
tering of otherwise free electrons, causing a phase shift 
5p = 7r/2 of their wave functions [cf. Eq. I|15f) ]. It is in- 
structive to make this duality transparent by identifying 
the effective low-energy Hamiltonian H c h arge that pro- 
duces the charge spectrum in Eq. 111|) . Introducing the 



the absence of any relevant operator [in the renormaliza- 
tion group (RG) sense] that could take the model to a 
new boundary fixed point - and hence change the gluing 
conditions in Eqs. i|12|) - we expect that the finite-size 
spectrum still has a Kac-Moody structure with Fermi liq- 
uid gluing conditions. Since the groundstate now has an 
(average) surplus charge n c ^ 1 compared to the case 
with an empty localized level, we conclude from Eq. @ 
that 



(13) 

I 

changed), will be shown to be self-consistent when we 
analyze the full two-channel Anderson model. 

We can corroborate our procedure by comparing 
our result in Eq. I|13|) with that obtained by Fuji- 
moto and collaborators^ Applying standard finite-size 
techniques^ to the exact Bethe Ansatz solution^! of the 
Anderson model these authors found that 



+ 0(±), (14) 

I 

charge currents 

J{x) =:^{x)i> a(T {x): (16) 

(where as before the normal ordering is taken with re- 
spect to the filled Fermi sea) , one immediately recognizes 
the charge part of lfTT|) as the spectrum of the Sugawara 
Hamiltonian^ 

Hcharge = \ J dx : J{x)J{x) ! - 1 J(0) . (17) 

The term ~ J(0) explicitly reveals the presence of an 
effective local scattering potential. By redefining the 
charge quantum numbers in the integer valence limit, 
Q — 1 — > Q, this potential is disguised as a renormal- 
ized boundary condition, encoded in the Kondo gluing 
conditions I|1UI) . 

With these preliminaries let us now go back to the full 
two-channel Anderson model in Eq. first considering 
the case when e — e s — e q <C /i — T. From the Bethe 
Ansatz solution— one finds that in this limit the physics 
is that of the overscreened two-channel (magnetic) Kondo 
models This is expected, since for this case the ion has 



(j(e-»j 2 -(^) s +iv c ) + (^%ti) + iv.) 
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integer valence, with an associated magnetic moment of 
spin 1/2. The finite-size spectrum of the two-channel 
Kondo model has been derived from that of free electrons 



(carrying spin and flavor) via Kac-Moody fusion*^ with 
the spin-1/2 conformal tower One finds, to 0(\/i), 



J 



E = En + - 



js(js + 1) 



+ N f 



(18) 



r 



Here Q G Z are U(l) charge quantum numbers, while 
j s (ff — 0,1/2,1 are quantum numbers for the level-2 
SU(2) spin (flavor) primary states, with N C) N s , and Nf 
labeling the corresponding descendant levels. The quan- 
tum numbers are restricted to appear in the combinations 



Q = mod 2, 



3 s 



j f = or 1, 
1 



Q = lmod2, j s =0 orl, j f = -, (19) 

which define the two-channel Kondo gluing conditions 19 
for the conformal towers. Since the charge quantum num- 
bers in Eqs. (|19fl are defined modulo 2, we can make a 
shift Q — > Q — 2 in (|18fl without affecting the gluing 
conditions. With this we recover the form of the confor- 
mal spectrum as derived by Fujimoto and Kawakami 34 
from the exact solutionrfS2& By mapping the nontrivial 
Z 2 part of the two-channel Kondo scattering onto that 
of a restricted solid-on-solid model 3 ^ coupled to the im- 
purity, these authors elegantly circumvented the "string 
solution" problem mentioned above. 

In exact analogy with the single-channel case detailed 
above, we can match the finite-size two-channel Kondo 
spectrum defined by by Eqs. I|18fl an d i|19|) to that for 
the n c — ► 1 limit of the two-channel Anderson model by 
performing a shift Q — > Q — 1 in Eqs. lfT%|) and ijPbil) . The 
gluing conditions are accordingly modified and now take 
the form 

1 



Q = mod 2, ] s = 0, or 1 
Q = 1 mod 2, j s = i, j f = or 1. (20) 



Jf = 5> 



The gluing conditions 120f) can formally be obtained by 
starting with the gluing conditions for free electrons (car- 
rying spin and flavor): 



Q = mod 2, 



Q = 1 mod 2, j s = i, j f = - 



0, jf = 0, 

l 

2' 



and then performing Kac-Moody fusion with the jf = 1/2 
flavor conformal tower: jf = — > 1/2; jf = 1/2 — > 
0,1; jf = 1 — * 1/2. Alternatively, we can obtain Eqs. 
(|20|l from Eqs. I|21|l by fusion with the j s — 1/2 spin 
conformal tower: j s = — > 1/2; j s = 1/2 — > 0,1; j s = 
1 — > 1/2, concurrent with fusion with the Q — 1 charge 
conformal tower: Q — > Q + l. In Sec. (III.B) we shall see 
how these two equivalent schemes can be given a simple 
physical interpretation, considering the particular struc- 
ture of the hybridization interaction in 10}. It is impor- 
tant to realize that in contrast to the single-channel case, 
the two-channel Kondo gluing conditions in l|19|) cannot 
be traded for a local scattering potential by the charge 
shift Q — > Q — 1. Although a local scattering potential is 
generated [cf. Eq. (|23(l below] the "charge-renormalized" 
gluing conditions 120|) are not those of free electrons. A 
two-channel Kondo system forms a non-Fermi liquid with 
properties very different from those of phase-shifted free 
electrons. 



(21) 



Let us now turn to the case of noninteger impurity va- 
lence n c ^ 1. As supplementary input we use the fact 
that the low-lying excited states of the model split into 
separate gapless charge, spin and flavor excitations for 
any value of e (and hence n c ), as revealed by the ex- 
act Bethe Ansatz solution^ Since there is no change of 
symmetry of the model as we tune e q and e s away from 
the (magnetic two-channel Kondo) limit e -C /i — T, one 
expects that a spectral U(l) charge, SU(2)2 spin, and 
SU(2)2 flavor Kac-Moody structure is retained through- 
out the full range of magnetic and quadrupolar energies. 
This implies that in the absence of any relevant operators 
(that could push the model to a fixed point of a differ- 
ent Kac-Moody structure) , the finite-size spectrum must 
take the form 



E = E a + - 



\{Q-n c f-\nl 



N, 



js (js + 1) 



N, 



(22) 
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with the values of the charge- spin- and flavor- quantum 
numbers constrained by Eqs. i|20|) . Note that the charge 
levels have again been shifted by a constant so that Q 
measures the number of conduction electrons added to 
or removed from the groundstate [i.e. Eq = E(Q = 0)]. 

Similar to the ordinary Anderson model, the charge 
spectrum in IL'l't implies the presence of an effective local 
scattering potential 

V*ff = ~J(0). (23) 

The charge current J has scaling dimension A = 1 , mak- 
ing V e ff in Eq. I|23|) into an exactly marginal bound- 
ary operator with the scaling field sampling the impu- 
rity charge valence n c . Provided that no relevant op- 
erators intervene under renormalization, it produces a 
line of stable fixed points parameterized by the value of 
n c or, equivalently, e. This is in accordance with the 
Bethe Ansatz solution^ which established a line of low- 
temperature fixed points analytically connected to the 
overscreened two-channel (magnetic) Kondo fixed point 
at e <C n — r. 

Given the finite-size spectrum defined by Eqs. I|2U[I and 
(|22[) we shall next identify the leading operator content 
of the scaling theory, following the "double-fusion" pre- 
scription outlined in Sec. II. In particular, we must check 
that no relevant operators are generated as we vary e and 
move away from the integer valence limit n c = 1. 

Performing a second Kac-Moody fusion- with the 
jf = 1/2 flavor conformal tower, the gluing conditions 
(|2T)j> change to 

Q = mod 2, j s — or 1, jf = or 1, 

= 1 mod 2, j s = \, if = \- ( 24 ) 

These new gluing conditions ensure that the boundary 
condition in the half-plane geometry is time independent 
(corresponding to having attached an Anderson impu- 
rity to each edge, x = and x = £, of the strip; cf. 
Sec. II). As we are now effectively considering a finite- 
size spectrum with two quantum impurities present, the 
charge quantum number gets renormalized twice, with 
Q — > Q — n c (x = 0) = Q' from the impurity at the x = 
edge, and Q 1 — > Q' — n c (x = £) = Q' + n c (x = 0) = Q 
from the x = I edge. Thus, the modifications of Q caused 
by n c at each edge of the strip cancel each other, and the 
double-fusion spectrum becomes independent of the im- 
purity valence. The reason for the cancellation is the 
same as in the single-channel case, where the sum rule 
(|15|l connects n c to the scattering phase shift caused by 
the impurity. Since the phase shifts at x = and x = £ 
have opposite signs (at one edge a left-moving electron is 
reflected, while at the other edge a right-moving electron 
gets reflected), one must formally assign charge valences 
of opposite signs to the auxiliary impurities attached to 
x = and x = £, respectively. We clearly expect the same 
situation to apply also in the two-channel case. By com- 
parison with Eqs. © and i|18|) we thus read off for the 



possible boundary scaling dimensions A = A c + A s + Aj : 

a c = ig 2 + A c , 

A s = \js{]s + 1) +N S , A f = lj f (jf +1)+N f 

(25) 

with Q, j s , and jf constrained by 124|) . and with 
N c , N S7 Nf positive integers. 

The gluing conditions in Eqs. I|24|) are identical to 
those for the spectrum of boundary scaling dimensions 
in the two-channel Kondo model. 19 Here we obtained 
them by double fusion with the jf = 1/2 flavor con- 
formal tower, while in the Kondo case they follow from 
double fusion with the spin-1/2 conformal tower. The 
equivalence of the two procedures is consistent with the 
Qmod2 invariance of Eqs. (|24|l . Double fusion with 
jf = 1/2 in the empty valence limit turns into double fu- 
sion with j s = 1/2 in the Kondo limit with a concurrent 
charge renormalization Q — > Q + 2. We shall elaborate 
on this point in the next section. 

The symmetries of the Hamiltonian in Eq. Q impose 
further restrictions on the allowed Kac-Moody charge, 
spin and flavor quantum numbers. Charge conservation 
implies that Q = in Eq. (|25l) and hence only the trivial 
identity conformal tower in the charge sector contributes 
to the spectrum of boundary scaling dimensions. More- 
over, since the charge valence n c does not appear in Eqs. 
(|25|1 . no scaling dimension can depend on it. No relevant 
operator can therefore appear as n c is tuned away from 
the magnetic two-channel Kondo limit. This guarantees 
that the finite-size spectrum defined by Eqs. I|24() and 
(|25|l remains valid throughout the full range of charge 
valence n c £ (0, 1), with a line of stable fixed points 
parametrized by n c . Recall that this line is produced 
by the marginal operator V e ff in Eq. I|23() . associated 
with the "single- fusion" spectrum, Eqs. (|2*U1) and J22J, 
which carries an explicit dependence on n c . V e ff here 
plays the role of a boundary changing operator^ For a 
given value of n c it determines the locus on the fixed 
line to which the theory flows under renormalization. It 
may be tempting to conclude that the result obtained 
from double fusion, Eqs. I|24|) and l|25|) . implies that the 
time-independent critical behavior is the same along the 
critical line, being insensitive to n c . This is correct with 
regard to critical exponents, but as we shall see below, 
amplitudes of various scaling operators pick up a depen- 
dence on n c , giving rise to very different physics as we 
move along the fixed line. 

The leading boundary operator contributed by the 
charge sector is the exactly marginal charge current 
J(0) =: ^ CT (0)^ Qff (0) :. It is formally identified as 
the first Kac-Moody descendant J(0) = J_il (c) (0) 
of the identity operator (which generates the Q — 
conformal tower). As we have seen, J(0) being 
marginal, it shows up manifestly in the low-energy 
spectral-generating Hamiltonian H charge in Eq. i|17fl . 
and is responsible for producing the line of fixed points. 
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The appearance of the charge current breaks particle-hole 
symmetry, since J(0) — * — J(0) under charge conjugation 
ip aC r(x) — > CcpEv^ipp (x). It is interesting to note that 
while the same symmetry is broken by the ionic pseu- 
doparticle term © in under /J — > e a ^f^, the pieces 
(0 and (@J of the Hamiltonian involving the conduction 
electron field %p arT (x) actually respect charge conjugation. 
This symmetry gets broken dynamically via the coupling 
to the pseudo-particles, allowing the charge current in 
Eq. (|2S1) to enter the staged 

The next-leading boundary operator from the charge 
sector is that of the energy-momentum tensor, 
T< c )(0) = : J(0)J(0):/4, appearing as the second Vi- 
rasoro descendant of the identity operator, T c (0) = 
L_ 2 l (c) (0)- Although it carries scaling dimension A = 2 
and hence is subleading to the charge current J(0), being 
a Virasoro descendant of the identity operator it has a 
nonvanishing expectation value and turns out to produce 
the same scaling in temperature as J(0). Although not 
central to the present problem, we shall briefly return to 
this question in Sec. IV. 

Focusing now on the spin and flavor sectors, conser- 
vation of total spin and flavor quantum numbers re- 
quires that all spin and flavor boundary operators must 
transform as singlets. The leading spin boundary op- 
erator with this property, O^(0) call it, has dimen- 
sion A s = 3/2 and is obtained by contracting the spin- 
1 field 4>^ s \0) (which generates the j s = 1 confor- 
mal tower) with the vector of SU(2)2 raising operators 
j« . W(o) = j(_ s ) . 0( S )(O). As expected, this is the 
same operator that drives the critical behavior in the two- 
channel (magnetic) Kondo problem^ The next-leading 
spin boundary operator is the energy momentum tensor 
T«(0) =:J (s) (0) • J (s) (0): /4 [= L_ 2 1«(G)], of dimen- 
sion A s = 2. 

The two leading flavor boundary operators are ob- 
tained in exact analogy with the spin case. In obvious 
notation: O^(0) = J U \ (0) of dimension A } = 3/2, 

andT</)(0) = : J (/) (0)- J (/) (0) : [= L_ 2 l (/) ( )] of dimen- 
sion Af — 2 . 

Boundary operators with integer scaling dimensions 
generate an analytic temperature dependence of the im- 
purity thermodynamics, 20 subleading to that coming 
from the two leading irrelevant operators O^ s \0) and 
0(^(0). Hence, in what follows we shall focus on these 
latter operators. In the case of the two-channel Kondo 
problem the flavor operator O^^(0) was argued to be ef- 
fectively suppressed, at least for the case when the (bare) 
Kondo coupling is sufficiently smalltA^ On dimensional 
grounds one expects that the spin scaling field X s [multi- 
plying 0^(0) in the scaling Hamiltonian] is 0(1/s/Tk), 
where Tk is the Kondo scale for the crossover from weak 
coupling (high-temperature phase) to strong renormal- 
ized coupling (low-temperature phase). The flavor oper- 
ator, on the other hand, does not see this scale since the 
infrared divergences in perturbation theory (which signal 
the appearance of a dynamically generated scale) occur 



only in the spin sector. The only remaining scale is that 
of the band width D (which plays the role of an ultravio- 
let cutoff), implying that the flavor scaling field A/ [mul- 
tiplying {f) (Q)] is 0{1/Vd). For a small (bare) Kondo 
coupling A, Tk ~ Dexp(— 1/A) <C D, and the critical 
behavior is effectively determined by (0) alone. 

As we will show, the picture changes dramatically for 
the two-channel Anderson model. The Bethe Ansatz 
solutionis shows that there are two dynamically gen- 
erated temperature scales T±(e) present in this model. 
They determine and parametrize the thermodynamic re- 
sponse of the system. Consider for example the quench- 
ing of the entropy as the temperature is lowered: in the 
magnetic moment regime (e <C fJ> — T) the two scales 
are widely separated, T+ <C T_. For temperatures in 
the range T + < T < T_ charge and flavor fluctuations 
are suppressed and the entropy shows a plateau corre- 
sponding to the formation of a local moment. The mo- 
ment undergoes frustrated screening when the temper- 
ature is lowered below T_ reaching zero point entropy 
S = fcsln-s/2, and n c = 1. As £ decreases T + and 
T_ approach each other and the range over which the 
magnetic moment exists decreases too. It disappears 
(i.e. T + = T-) when e = \i. At this point no mag- 
netic moment forms. However the system still flows to a 
frustrated two-channel Kondo fixed point, with entropy 
S = ks In V2, but with n c = 1/2. Continuing along the 
critical line, the two scales trade places, and eventually, 
at the quadrupolar critical endpoint (e [i — T), one 
finds that T + 3> T_ . Since the irrelevant spin (flavor) op- 
erator is expected to dominate the critical behavior at the 
magnetic (quadrupolar) fixed point, we are led to conjec- 
ture that the corresponding scaling field is parametrized 
precisely by 1/ \/T+ (1/ \fTZ) . With this scenario played 
out, the relative importance of the two boundary oper- 
ators changes continuously as one moves along the fixed 
line, with each operator ruling at its respective critical 
end point. 

In the next section we shall prove our conjecture by 
computing the impurity specific heat and susceptibility 
following from the scaling Hamiltonian and comparing it 
to the low temperature thermodynamics from the exact 
Bethe Ansatz solution.— We shall also use this compari- 
son to "fit" the parametrization of the scaling fields. This 
fixes the critical theory completely. 



B. Application: Fermi Edge Singularities 

Before taking on this task, however, we shall exploit 
the BCFT scheme developed above to determine the scal- 
ing dimensions of the pseudoparticles which enter the 
Hamiltonian in and The analogous problem 

for the ordinary Anderson model has been extensively 
studied22ilMi, in part because of its close connection to 
the X-ray absorption problem. 

One is here interested in the situation where X-ray 
absorption knocks an electron from a filled inner shell 
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of an ion in a metal into the conduction band. As an 
effect the conduction electrons experience a transient lo- 
cal potential at the ion which lost the core electron, and 
this typically produces a singularity in the X-ray spec- 
trum ~ (oj — lof)~P close to the Fermi level threshold 
ujp. Nozieres and de Dominicis^ using a simple model, 
showed that the exponent depends only on the phase 
shift which describes the scattering of conduction elec- 
trons from the core hole created by the X-ray. Fermi 
edge singularities linked to local dynamic perturbations 
are fairly generic^ and it is therefore interesting to ex- 
plore how they may appear in a slightly more complex 
situation where a localized "deep hole" type perturba- 
tion involves additional degrees of freedom. Indeed, the 
two-channel Anderson model in Eq. Q provides an ideal 
setting for this as it accommodates a local quantum im- 
purity carrying both spin and channel (quadrupolar) de- 
grees of freedom. 

To set the stage, let us make a Gedankenexperiment 
and imagine that we suddenly replace a thorium ion in 
the host metal (say, Uo.gTho.iBeis) by a uranium ion, 
the sole effect being to introduce a localized quadrupolar 
or magnetic moment. To describe the response of the 
host to this perturbation, we consider the pseudoparticle 
propagators 



G a = (0\e iMt f a e- iHt fl\0) 



(0\e^%e- lH %\0) 



(26) 
(27) 



(with no summation over a, a). At time t < the ground 
state 1 0) of the metal is that of free conduction electrons. 
At t — a pseudoparticle (creating the localized mo- 
ment) is inserted into the metal, and the system evolves 



in time governed by the two-channel Anderson Hamilto- 
nian in Eq. (JIJ . At some later time t the pseudo-particle 
is removed, and one measures the overlap between the 
state obtained with the initial unperturbed groundstate. 
It is important to realize that by inserting (removing) 
a pseudoparticle operator, one simultaneously turns on 
(off) the impurity terms in Eqs. J3J and and hence 
changes the dynamics of the conduction electrons. As we 
reviewed in Sec. (II), at low temperatures and close to a 
critical point, this corresponds to a change of the bound- 
ary condition which emulates the presence (absence) of 
the impurity. Hence, the pseudoparticle operators f a and 
b a are boundary changing operators^ and we can use 
the BCFT machinery to calculate their propagators in 
Eqs. i|26l) and (|27|l . Our discussion closely follows that 
in Ref^H 

We start by considering the free-electron theory de- 
fined on the complex upper half-plane C + = {z — 
t + ix;x > 0} with a trivial (Fermi liquid) boundary 
condition, call it "A" , imposed at the real axis x = 0. 
We denote by | A; 0) the ground state for this configura- 
tion. By injecting a pseudoparticle at time r = 0, the 
boundary condition for r > changes to "£?", here la- 
beling the boundary condition which corresponds to the 
nontrivial gluing condition in (|20|l for the spectrum of the 
two-channel Anderson model. By mapping the half-plane 
to a strip {w = u + iv; < £} via the conformal trans- 
formation w — (£/2 / n)\nz 1 the boundary of the strip at 
v = £ also turns into type B after insertion of the pseu- 
doparticle operator. Under the same transformation the 
pseudoparticle propagator in C + , (A; | Oi(ri)0\(T2) \ 
A: 0) = (t\ — T2)~ 2xi (with xi the dimension of Oi,i = f 
or b) transforms as 



(n - r 2 y 



21 . TT 

— sinn — [ui —uo, 
7T 2£ y 



-2xi 



= J2\ <^;0|Oi(0) | AB;n) | 2 e 
I 



-(B^ B -E^ A )(«2-«i) 



(28) 



where on the right-hand side of Eq. (|28() we have inserted 
a complete set of states | AB; n),n = 0, 1, ... of energies 
E^ B , defined on the strip with boundary condition A (B) 
at u = (u = £). In the limit «2 - «i > I'. 



- S mh-( Ml 



u 2, 



^—[lTXi(u2-Ul)/g] 



It follows, by comparison with Eq. I|28l) . that 



7T 



AB 



E, 



AA\ 



(29) 



(30) 



where Ui is the lowest energy state with a nonzero matrix 
element in the sum in Eq. I|28|l . 

The pseudoparticle spectrum Ii(u>) close to the Fermi 
level is obtained by Fourier transforming the correspond- 
ing propagator in Eq. (|26|) or (|27|l . and it follows that 



■ LOF 



11-2X4 ' 



i = f,b, 



(31) 



with a singularity when Xi < 1/2. 

Let us consider first the slave boson b a , carrying quan- 
tum numbers Q = 0, j s = 0,jf = 1/2. Its scaling di- 
mension Xf, is given by Eq. (|30f) where E^ B is obtained 
from Eq. (2U by putting Q = 0,j s = 0, j f = 1/2, and 
with Eq A = (with proper normalization of energies). 
To be able to compare energies corresponding to different 
boundary conditions it is important to keep the overall 
energy normalization fixed, and hence we must remove 
the normalization constant ~iml/8£ in Eq. J^J before 
reading off the answer. (Recall that this normalization 
constant is expressly designed to remove the Q = con- 
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tribution from the charge sector to the ground state en- 
ergy.) We thus obtain 

Xb = h + \ n * ■ (32) 

Similarly, we obtain for the pseudo fermion (Q — 1, j s = 
l/2,i/ = 0): 

It is interesting to compare Eqs. (EU and with the 
corresponding scaling dimensi ons in t he ordinary Ander- 
son model, obtained in Refs. f39l41t x' b — n^,/2,x'f = 
1/2 — n c /2 + n^/A. We conclude that by opening an 
additional channel for the electrons, the singularity in 
the slave boson spectrum in Eq. (|31|) (corresponding to 
X-ray photoemission) gets softer, whereas for the pseud- 
ofermion spectrum (corresponding to X-ray absorption) 
the opposite is true. One can show that this trend is 
systematic, and gets more pronounced as the number of 
channels increases^ We also note that the values of the 
pseudoparticle scaling dimensions as obtained by various 
approximation schemes, e.g. the "non-crossing approx- 
imation" (NCA) and large-A^ calculations, 46 are not in 
agreement with the exact results presented here. 

Before closing this section, let us return to the 
Gedankenexperiment where we inserted an empty (de- 
generate) impurity level by adding a slave boson to the 
filled Fermi sea. The slave boson carries only flavor 
jf = 1/2, and it follows that the allowed values of 
the quantum numbers of the combined system are ob- 
tained by coupling jf = 1/2 to the flavor quantum num- 
bers j'f of the conduction electrons (leaving charge and 
spin untouched). Since jf and j'f label conformal tow- 
ers, this coupling is precisely governed by conformal fu- 
sion with Jf = 1/2 : j' f -| j' f - 1/2 |,| j' f - 1/2 | 
+1, ... min(jj + 1/2, 3/2 — j'f). By letting the system re- 
lax, a non-zero charge valence n c may then build up at 
the impurity site (as determined by the particular value 
of e — s s — e q ). Formally, this process is captured by 
the charge renormalization Q — ► Q — n c in the finite size- 
spectrum <|22[) . taking place concurrently with the fusion 
in flavor space which produces the nontrivial gluing con- 
ditions (|20|l . Running the Gedankenexperiment with a 
pseudo-fermion instead (carrying charge and spin, but no 
flavor) leads to a scenario dual to the one above, where 
the coupling is now governed by fusion with j s — 1/2 in 
the spin sector and Q = 1 in the charge sector [cf. the dis- 
cussion after Eq. (|21|l ] . As the system relaxes an average 
charge n c goes to the impurity site, again leading to the 
renormalization Q — ► Q — n c in the spectrum (|22p. We 
conclude that our key results, Eqs. i(2*0"|) and 1(2*2^1 - ob- 
tained in Sec (II. A) via a formal analysis - can be given a 
natural interpretation by considering the pseudoparticle 
structure of the hybridization interaction in Eq. Q). 



IV. LOW-TEMPERATURE 
THERMODYNAMICS 

A. Zero- Temperature Entropy 

The Bethe Ansatz solution 17 of the model reveals 
that the line of fixed points is characterized by a zero- 
temperature entropy Si mp = k B In \[2 typical of the two- 
channel Kondo fixed point. As we shall show next, this 
property emerges naturally when feeding in our result 
from Sec. (II. A) into the general BCFT formalism devel- 
oped by CardyiSi 

Recall that the key ingredient in this formalism is the 
modular invariancel* of a conformal theory. Applied to 
a model defined on a cylinder of circumference (3 = 1/T 
and length £, and with conformally invariant boundary 
conditions A and B imposed at the open ends, this means 
that its partition function Z AB is invariant under ex- 
change of space and (Euclidean) time variables. In our 
case, B is the boundary condition obtained from A by 
attaching an Anderson impurity to one of the edges of 
the cylinder (with A a trivial Fermi liquid boundary con- 
dition, to be defined below). Now, let H* be the critical 
Hamiltonian that corresponds to the spectrum given by 
Eqs. CD and (23. Then H AB = (£/n)H* generates 
time translations around the cylinder specified above. 
The partition function Z AB can hence be written 

Z AB = Tr e-^ B = ]T n a ABXa {e-* ' 1 ) . (34) 

a 

Here \a is a character of the conformal U(l) x SU(2)2X 
SU(2)2 algebra, with a — (Q,j s ,jf) labeling a product 
of charge, spin, and flavor conformal towers [constrained 
by Eqs. (|21l) ]. and with n a AB counting its degeneracy 
within the spectrum of H A b- By a modular transforma- 
tion (3 <-> I the role of space and time is exchanged, and 
the partition function now gets expressed as 

Z AB = {A\e- tH »\B), (35) 

where Hp is the transformed Hamiltonian generating 
translations along the cylinder, in the space direction, 
with | A) and | B) boundary states^ implicitly defined by 
Eq. I|35|) . By inserting a complete set of Ishibashi states 

I a ) = J2 m I a ' m ) L ® I a ' m )-R m t° Eq. (with m run- 

ning over all states in an L/R product of charge, spin, 
and flavor conformal towers labelled by the index a, one 
obtains^ 

Z AB = ^2(A\a) (a | B) Xa (e-^ e ^) . (36) 

a 

Equating Eqs. I|34|l and l|3t)|l one can immediately read 
off Cardy's equation*^ 

J2 n A B Sa b = (A\a)(a\B), (37) 

b 

where S a b is the modular S matrix defined by 

Xa(e-^ e ) = Y, S ^ b (e' M/f3 ), (38) 

b 
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with a and 6 indexing products of charge, spin, and flavor 
conformal towers. It follows that the nontrivial boundary 
state | B) (which emulates the presence of the impurity) 
is related to | A) via the identity^ 



(a\B)=(a\A) 



Sad 
SaO 



(39) 



with labeling the product of identity conformal tow- 
ers and d the product of conformal towers with which 
the Kac-Moody fusion of the finite-size spectrum is per- 
formed. As we found in Sec. (II. A), for the two-channel 
Anderson model, d = (Q = 0,j s =0,jf = l/2). Factoring 
the modular S matrices in charge, spin, and flavor, with 
a = (Q>3s,jf) constrained by Eqs. if^TJl but otherwise 
arbitrary, one obtains 



Sad 
SaO 



f 



QC OS C 

D Q0 D j s D j f l/2 



S 



QC QS Q f 



f 

3 1 1/2 

S f 
'j/0 



(40) 



The modular S matrix for the SU(2)2 flavor symmetry is 
given bji 2 * 



Q 1 «i 

iJij-q' - — ^ Sill 

V2 



n(2j f + l)(2j' f + l) 



and it follows from Eqs. and (gUJ) that 



(a\B) = (a\A) 



sin[f(2j/ 
sm[f(2i/ 



ill 
1)] 



(41) 



(42) 



We now have the tools for deriving the impurity en- 
tropy. Following the analogous analysis of the multichan- 
nel Kondo model in Ref . 1471 we take the limit £//3 — * oo in 
Eq. (|36[1 (i.e. we take /3 — > oo after the large- volume limit 
I — > oo). As a result, only the character of the ground- 
state, Xa=(o,o,o)(e~ 47r£/ ' 3 ) = e^ 6 ^, contributes to Z AB 
(with c the conformal anomaly of the theory) . One thus 
obtains 



Z AB ^e^ c ^(A\0)(0\B), t/(3- 



oo, 



(43) 



with | 0) =| Q = 0,j s = 0,jf = 0). The free energy is 
thus 



Fai 



-irck B T 2 £/6 - k B Thi((A | 0)(0 | B)). (44) 



The second term in Eq. (|44|l is independent of the size 
of the system, and we therefore identify 



Si; 



k B \n((A\0){0\B)) 



(45) 



as the impurity contribution to the zero-temperature en- 



tropy. With no impurity present, A = B and SV, 



0, 



and it follows that | A) = | 0) [consistent with the Fermi 
liquid gluing condition in Eq. 121|) ]. Thus, from Eqs. 142|) 
and H45fl we finally get 



Simp = k B In a/2 , 



(46) 



in agreement with the Bethe Ansatz result in Ref. ITtI 

A few remarks may here be appropriate. First, we note 
that the effective renormalization of the charge sector, 
Q — > Q — n c [cf. Eq. 1)22(1 ]. does not influence the impu- 
rity entropy. At first glance, this may seem obvious: the 
degeneracies of the impurity levels — and hence the im- 
purity entropy — is a priori not expected to be affected 
by the potential scattering in (|23|l caused by the impurity 
valence n c . The picture becomes less trivial when one re- 
alizes that the impurity valence in fact controls how the 
spin and flavor degrees of freedom get screened by the 
electrons. As found in the Bethe Ansatz solution^, for 
nonzero e = e s — e q , the quenching of the entropy oc- 
curs in two stages as the temperature is lowered. For 
sufficiently large values of | e — /i | , the free impurity en- 
tropy k B In 4 is first reduced to k B In 2 at an intermediate 
temperature scale, suggestive of a single remnant mag- 
netic or quadrupolar moment (depending on the sign of 
e or, equivalently, the sign of n c — 1/2). As the temper- 
ature is lowered further, this moment gets overscreened 
by the conduction electrons, as signaled by the resid- 
ual entropy k B In \/2. Remarkably, as found in Ref. IT^ 
and verified here within the BCFT formalism, the same 
residual impurity entropy k B ln^/2 appears also in the 
mixed- valence regime | s — fi | < T (where n c « 1/2), for 
which there is no single localized moment present at any 
temperature scale. 



B. Impurity Specific Heat 

We now turn to the calculation of the impurity spe- 
cific heat. Again, the analysis closely parallels that for 
the two-channel Kondo model,- 9 Some elements are new, 
however, and we here try to provide a self-contained 
treatment. 

First, we need to write down the scaling Hamiltonian 
H SC aiing that governs the critical behavior close to the 
line of boundary fixed points. To leading order it is ob- 
tained by adding the dominant boundary operators found 
in Sec (III) to the bulk critical Hamiltonian H* [repre- 
senting H bulk in Eq. JJJ]: 

H scalmg = H* + A c J(0) + \ s O^{0) + \fOW(0) 

+ sublcading terms , (47) 

with A C , S) / the corresponding conjugate scaling fields. As 
reflected in the finite-size spectrum (I22|l . H* splinters 
into dynamically independent charge, spin, and flavor 
pieces, implying, in particular, that all correlation func- 
tions decompose into products of independent charge, 
spin and flavor factors. On the critical line the impu- 
rity terms Hi on and Hhybr in Eqs. and Q together 
masquerade as a boundary condition on H* , coded in 
the nontrivial gluing conditions in Eqs. (|20|) . In addition 
Hi on and Hhybr also give rise to the exactly marginal 
charge current term A c J(0) in Eq. (|47|l . This term keeps 
track of the charge valence acquired by the impurity. As 
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we have already noticed, this is different from the Kondo 
problem where the local spin exchange interaction at crit- 
ically gets completely disguised as a conformal boundary 
condition^ The reason for the difference is that away 
from the (integer valence) Kondo limit, the effective lo- 
cal scattering potential in Eq. (|23|) cannot be removed 
by redefining the charge quantum numbers Q — > Q — 2 
in Eqs. JSHJl and This in turn reflects the fact that 

for non-integer valence the dynamics manifestly breaks 
particle-hole symmetry, whereas for integer valence there 
is an equivalent description^ in terms of the two-channel 
Kondo model where this symmetry is restored^. Off the 
critical line the impurity terms in Eqs. © and (J2J ap- 
pear in the guise of irrelevant boundary terms, of which 



A S 0( S )(O) and A/OW)(0) in Eq. are the leading ones. 

Recall from Sec. (II) that a boundary critical the- 
ory by default is defined in the complex half-plane 
C + = {Imz > 0}, with the boundary condition im- 
posed on the real axis Imz = 0. As we have seen, 
by continuing the fields to the lower half-plane we ob- 
tain a chiral representation of this theory, now defined 
in the full complex plane C. Via the transformation 
w = t + ix = {(3/ir) arctan(z), C can be mapped confor- 
mally onto an infinite cylinder, T c call it, of circumference 
0: T c = {(t,x);-(3/2 ^ t < (3/2; -oo s$ x < oo}. It is 
convenient to use T c as a "finite-temperature geometry" 
by treating (3 as an inverse temperature. The partition 
function can then be written as 



r0/ 2 

e -0F(0,{K}) = e -0F(0,o) ( exp ( / oIt[\ c J(t, 0) + \ s 6 {s) (r, 0) + \ f 6 {f) (r, 0) ] + subleading terms ) ) T , (48) 

J -0/2 

where F((3,{X U }) — Fb u ik((3,0) + fimp(0,{X v }) is the sum of bulk and impurity free energies, with {A^} = 
(A c , A s , A/, ...) the collection of boundary scaling fields. We have here passed to a Lagrangian formalism, with the 
"tilde" and ( )t referring to the finite-T geometry T c . By a linked cluster expansion we obtain from Eq. (|48|l the 
impurity contribution to the free energy: 

fbnp(P,M) = -Za / / dTidT 2 [X 2 c (J(n,0)J(T 2 ,0)) T ,c+ V A 2 , (O u ( n , 0)O»(t 2 , 0)) t ,c ] + subleading terms , 

, (49) 



with ( )t, c denoting a connected two-point function in 
T c and where, for a fixed temperature 1/(3, the constant 
term fi mp ((3, 0) has been subtracted. Note that the linear 
terms in the expansion are absent since all three bound- 
ary operators are Virasoro primary 24 and hence have van- 
ishing expectation values. Since two-point functions con- 
taining operators from different sectors (charge, spin or 
flavor) decompose into products of one-point functions, 
these also vanish. Similarly, connected two-point func- 
tions reduce to ordinary two-point functions. Also note 
that there is no term in Eq. (|49|l of zeroth order in the 
scaling fields. This is a generic feature of quantum im- 
purity problems, which is best understood by looking at 
the finite-size scaling forrn^ for the impurity free energy: 



Thus, to see the effect of the impurity, the leading irrele- 
vant boundary operators have to be added explicitly, and 
these operators then produce the dominant terms in the 
scaling in temperature. 

The correlation function of chiral Virasoro primary op- 
erators in C takes the familiar form 



(0W(*i)0«(z a )> 



A, 



(zi - z 2 ) 2A * 



(52) 



(50) 



with Yi mp a universal scaling function and Ei mp nonuni- 
versal. The exponents y v are RG eigenvalues, connected 
to the dimensions A„ of the irrelevant boundary opera- 
tors by y v = 1 — A^. It follows immediately from Eq. H5()(l 
that at the fixed point 



Ci m p({A„ — 0}) — 



-T 



B 2 f- 

u J imp 

dT 2 



0. 



Here Aj is the scaling dimension of 0^ l \ and Ai is a 
normalization constant. Mapping the real axis Im z = 
into T c , this transforms into 



(0M(n,o)0«(7 a ,o)> T = - . * 

|fsm(-g(ri -7>!))| 



2A; 



(51) and it follows that Eq. (|49|l can be rewritten as 



(53) 



13 



fimp{P, {Ay}) — 



2(3 



du XzA, 



tans 



\ + { X 2 s A s + X 2 A f )^- 



subleading terms . 



(54) 



r 



with k a UV cutoff. This form is obtained by replacing 
the double integral in Eq. (|49|l by a single integral over 
t = t\ — T2, exploiting that the two-point functions in 
Eq. 153|) are periodic in r (with period 0). Introducing 
the auxiliary variable u = tan(7rr//3) and specifying the 
scaling dimensions A c = 1, A s = Af = 3/2 in Eq. (I5S 



one arrives at the expression in Eq. I|54|) after inserting 
the short-time cutoff tq = k(3/tt. The integral can be 
evaluated straightforwardly 19,49 and in the limit of low 
temperatures, i.e. for small k, an expansion of tan(«;) 
yields for the impurity specific heat: 



Cimp — 



-T 



d 2 f- 
dT 2 



= {X 2 A S + X 2 f 4/)i 2 Tln(— ) + X 2 c A c tt 2 t T 



corrections from subleading boundary operators. (55) 



r 



Before proceeding, let us comment on the correction 
terms in Eq. (|55|l . The next-leading boundary opera- 
tors are the energy-momentum tensors T^ c ' s '^ identified 
in Sec. (III. A). Passing to the finite-T geometry r c , their 
contribution to the partition function takes the form 



exp( £ / drX[f^(r 7 0))) T , 

,„ - * J -0/2 



i=c,a,f ' 



with A [ conjugate scaling fields. It follows, again by a 
linked cluster expansion, that they contribute 



Z P ■ rJ-B/2 
t—C,S,J 11 



drA' J (f«(r,0)) T 
+ subleading terms (56) 



to the impurity free energy. Note that in contrast to 
the contribution from the leading boundary operators in 
Eq. H49fl . fifimp is linear in the scaling fields A^. This 
reflects the fact that the energy-momentum tensors are 
Virasoro descendants of the identity operator, with non- 
vanishing one-point functions^ 



4 = 1, d a = df = ~. (57) 



Inserting Eq. l(57|l into (|5*dT) and integrating, it follows 
that the contribution to the impurity specific heat is 



3A; 
2 



3A' 



(58) 



This is the leading correction in Eq. (|55|l that is linear in 
the scaling fields: higher-order Virasoro descendants of 
the identity operator generate higher powers in tempera- 
ture. Similarly, higher-order Kac-Moody descendants of 
the cf>^ and <f>^ operators produce subleading correc- 
tions to the impurity specific heat that are quadratic in 
the scaling fields. 

Returning to Eq. (|55(l it is important to notice that the 
leading logarithmic terms are independent of the cutoff 
procedure. We can make this explicit by introducing the 
temperature scales T s and Tf that characterize the spin 
and flavor dynamics, and absorb the arbitrary cutoff r 
in the subleading terms: 



C im p = X 2 s A s n 2 TH^) + X 2 A f n 2 TH^) 



X 2 c A c t + X 2 S A S ln( 



t T s 



XjAfH—r) ) it T 



T Tf' 



subleading terms . (59) 



It remains to determine the parameters X s j,A s j and T s j that enter the leading terms in Eq. (|59|l . Start- 
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ing with the normalization constants of the two-point 
functions [cf. Eq. I|52[l ]. these are defined bj*24 A s j = 

(j)( s 'f)} m Dropping the indices (s, /), 

4>\ j-i ■ 0) = (r i i </>") = 



-1 

we have that 



jLY ; -0 (>,/) |j-i 



(*,/) 
j 



[Ji,j- 



where 



{(j> a | [Jf, J b _ x \ + JtxJx | 4> ) = 
the last identity follows from | (f> b ) being a Kac-Moody 
primary state. Using the SU(2)2 Kac-Moody algebra 
[Jf, J^x] = ie abc J§ + S ab (common for spin and flavor 
sectors), together with the fact that <j> transforms as a 
spin-1 object under global SU(2) (generated by Jo), one 
straightforwardly arrives at the result A s j — 9. 

Turning to the scaling fields A s j, on dimensional 
grounds these must take the form 



-W - 



B 



(60) 



with B s j dimcnsionless constants. The scaling fields and 
temperature scales T s and Tj (which also enter the loga- 
rithms explicitly in Eq. (jSHJl) will be extracted from the 
numerical solution of the Thermodynamic Bethe Ansatz 
(TBA). 17 But first we need to discuss the effect of applied 
external fields. 



must be added to the scaling Hamiltonian in Eq. I|47|l . 
Here J| = (1/2) '■^\ l J jZ u,v'^ ) av '■ i s the z component of the 
electron spin current that couples to the field. Passing 
to the finite-T geometry T c , this adds a term 



B 



dx 



<1tJ z s (t,x) 



to the impurity effective action in Eq. (|48|l . By a linked 
cluster expansion to 0(B 2 ), the impurity magnetic sus- 
ceptibility 



C. Magnetic Susceptibility 



d 2 f- 



dB 2 



5=0 



(61) 



In the presence of a magnetic field B, applied in the z 
direction, the term 



dxJ z (x) 



gets expressed as 



Ximp 



i=c,s,/ 



dxidx 2 



l=C,S,f 



d Tl dT 2 dT 3 {JlH6f) T ,c 



dx\dx 2 



dn...dT 4 {JlJ^dfd^T.c + subleading terms , (62) 



where = J, J z = J*,Jf = J z (r^x l ),of = 

O®(Tk,0) and = J { '[ 1] ■ (f> {sJ) . Note that there is 

no term to zeroth order in Ai, for the same reason that 
zeroth order terms were absent in the expansion of the 
impurity specific heat. Terms linear in A c and A/ are 
easily seen to vanish in 162|) due to charge-, spin-, and 



flavor decomposition of correlation functions: (O^) and 
(O^) factor out, and with these operators being Vira- 
soro primaries, the expectation values vanish. Writing 
out the remaining connected three-point function in the 
spin sector (going back to the complex plane for nota- 
tional simplicity), we have 



J 



(■MOW) 



(J z J z 0^)~(J?J z )(0i s) ) 



(63) 
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with J z = J'(z k ) and 0< fl) = 0^(z k ). A gain, all terms in Eq. (|63() containing one-point functions are zero since J z 
and 0( s ) are Virasoro primaries. Then, inserting the operator product expansions 50 

J z ( Zl )0 {s) {z 2 ) = j^ z (z 2 ) + regular terms in (zi - z 2 ), (64) 

(zi - z 2 y 

J z (zi)(f>^ z (z 2 ) — regular terms in (z\ — z 2 ), (65) 



into (Jf J 2 3 ), this three-point function collapses^ to 
a sum of vanishing one-point functions in the neighbor- 
hood of z\ ~ Z2 ~ Z3. Since the three-point functions in 
Eq. E3, (Jf JfO£ s) ) c and (Jf JfO^), both vanish for 
infinite separations of the fields, analyticity then implies 



I 

that ( Jf Jf03 )c = for any 21,22,2:3. One concludes 
that there is no contribution to Ximp linear in A s . 

Turning to the terms in Eq. i|62[l quadratic in the scal- 
ing fields, we have (again going back to the complex 
plane) 



(Jf J z 2 ofof) c = {J z x JiOfof) - {J{Jl){Ofof) - {J z x Of){J z 2 Of) - (J z O^)(J z oi l) )- 

— terms containing vanishing one-point functions . (66) 



As for the case above and as expected on physical 
grounds, only terms containing the spin boundary oper- 
ator 0^ s ' can potentially contribute in Eq. (|66|l . Inject- 
ing Eq. into {J[0 ( 3 s) ){J$0i s) ) and (Jf ( 4 s) )(J z O { 3 s) ), 
these terms are seen to vanish, and one is left with 

( Jf JiOPoP) (Jf J|) {0 { 3 S) { : ] ) . (67) 



The four-point function in Eq. (|67|l is easily calculated, 
again by exploiting the operator product expansions in 
Eqs. (|64|l and l|65|l . Transforming to the finite-T geome- 
try T c , one fmds»£i 



( r {t 1iXi )j z (t 2 , z 2 )0« (73, 0)6^ (7-4, 0)) 



T,c 



= 9 



h(Ti-T 3 +ixi)h(T 2 -T i +ix 2 ) + h(T 2 -T 3 +ix 2 )h(Ti-T A - 



^ sin ^(r 3 -r 4 ) r 1 , (68) 



where h{z) = [(/3/tt) sin(7rz//3)]~ 2 . Inserting Eq. (Sgjl 
into (|62() and carrying out the integration, one finally 
obtains, to O(A^) [with all other terms in Eq. (|62|l van- 
ishing to this order] : 

Xim P = l8\ 2 s Hj^) + 0(X 3 s ), (69) 

where we have subtracted a temperature-independent 
nonuniversal constant (what in practice fixes the scale 
T s ). Using the Bethe Ansatz result for T s and A s that we 
give in the next section, we finally obtain a complete de- 
termination of the asymptotic low-temperature behavior 
of the susceptibility. 

Notice that for the case of the response to an external 
field in the flavor degrees of freedom the calculations pro- 
ceed in a completely analogous way and the final result 



is the same with the obvious replacements T s — > Tf and 
A s -> A/. 



D. Numerical Fit 

We can use the numerical solution of the TBA equa- 
tions in Ref. Il7l to find the values of the conjugate scaling 
fields and temperature scales that parametrize the line of 
boundary conformal fixed points. This is, to our knowl- 
edge, the first time that this kind of fit is being explicitly 
carried out. 

We start from the following asymptotic expression for 
the free energy from which the above results for the en- 
tropy, specific heat, and susceptibilities follow immedi- 
ately: 
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F mp [{t, h sJ « t}] = i? mp - T + £ T 2 In ^ ( I 



x=sj 



T 

- 1 - r. 



7T 2 T 2 



(70) 



with = const, S^, = fc B In \/2, and /i S) / the external field acting on the spin or flavor degree of freedom 

respectively. In order to carry out the fit, we define the quantities 

2 _ . .„ _ „_ , n. T 



h 2 r 



(F imp [{T, h x « T, fee = 0}] - F imp [{T, h sJ = 0}]) = 18A X 2 In + O ({T 2 , fe^} , ) (71) 



r 



which coincide with minus the field susceptibilities in 
their low-T behavior and deviate from them as the tem- 
perature increases. We have here used the notation 
x = s (/) if x = f (s). The advantage of these quantities 
over the usual susceptibilities is that they incorporate the 
knowledge that we have about the form that the leading 
free energy terms take and they are numerically simpler 
to compute. 

We fitted the data from the numerical solution of the 
TBA equations in the appropriate temperature range us- 
ing the form xf mp ~ a-x^T + b x and extracted the de- 



sired parameters (as functions of e): X x = 




lnT x = -|^. We also computed B x = K^TT^. In Fig.Q 
we show the values obtained for A/ as a function of the 
doublets energy difference e (we have taken fj, = 0). 




FIG. 1: The dots are the Thermodynamic Bethe Ansatz nu- 
merical results for the flavor scaling field as a function of e 
(both in units of F) and the line is our exponential fit. 



In the same figure we show a fit using the function 




FIG. 2: The dots are the Thermodynamic Bethe Ansatz nu- 
merical results for the flavor temperature scale as a function 
of e (both in units of Y) and the lines are two different fits. 



also show in the same figure the results of fitting In Tf = 
a + b e which gives the values: a = 0.033 ± 0.012 and 
b = -1.581 ± 0.003. The fit was done for e > and 
the correlation coefficient was \R\ = 0.9999501. We next 
carried out the fit of a more complicated function in order 
to capture the behavior for e < as well. We fit using 
lnTy = a + b e + c e de where a and b are taken to be 
those determined above. The result of this fit is also 
shown in the same figure, and the values obtained for the 
parameters were c = 0.043 and d = —2.24. The fitting 
was done this time in the full range of e, with correlation 
coefficient \R\ = 0.998274. 

In the e > region, we combine the results of the two 
fits to obtain 

B f = Xfy/Tf w 0.1666 w 1/6. 



X f = X e 16 that gives the values: A = 0.1640 ± 0.0004 From the expression T t = that was derived in 

the Bethe Ansatz analysis^ we expect that b = —it/2 
in the region where e > T. This value is in fair agree- 
ment with the result of our fit carried out over a larger 
region (e > 0). The agreement not only validates our re- 
sults, but also serves as an indicator of the magnitude of 
the total errors of the numerical procedure (always nec- 
essarily larger than the statistical errors of the fit alone) . 



and 7 = 0.7838 ± 0.0004. The fit was done in the in- 
terval [—2,8] and the correlation coefficient was \R\ = 
0.9999935. The errors indicated are statistical ones and 
do not take into account other sources of error; the total 
errors are somewhat larger (as we discuss below). 

In Fig. |21 we show the values obtained this time for 
Tf as a function of the doublets energy difference. We 
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The reader should notice that the parameter a cannot 
be compared in a similar way since it is not a universal 
quantity and we did not adopt the same prescriptions in 
defining 7] and Tf (the prescription adopted in the latter 
case was that there are no constant terms in the defini- 
tion of xfmp an d that the subleading terms are regular, 
i.e. vanishing, functions of the temperature). 

Notice that for negative e the scale Tf grows exponen- 
tially instead of the linear growth of the scale associated 
with the quenching of charge fluctuations. This is no 
contradiction since the two scales are not expected to be 
equal on physical grounds. In other words, there is no 
reason to expect that the leading irrelevant operator in 
the flavor sector will govern the charge fluctuations when 
e < 0. One would rather expect that many operators in- 
cluding those in the charge sector of the theory would be 
involved, and in any case the quenching of charge fluc- 
tuations takes place at temperatures far removed from 
the domain of validity of the scaling Hamiltonian derived 
from BCFT. The exponential growth of Tf that we find 
is indicative of a very rapid suppression of the 'impor- 
tance' of this operator relative to its counterpart in the 
spin sector of the theory as e grows negative. 

Finally, we also studied the effects of applying a mag- 
netic field and confirmed numerically that 



X s (e - fi) = \f(fi- e) 
T s (e - n) =T f {^~ e) 



With this observation the parameters of the asymptotic 
form for the free energy are now fully determined for any 
value of e. 

A note on some of the details of the numerical calcu- 
lation is in order. The results for the free energy in the 
presence of an applied field were obtained using a con- 
stant field h s ,f — 10~ 10 r, which remained always at least 
two orders of magnitude smaller than the lower limit of 
the temperature range used in the fitting of Ximp- This 
was to ensure that we were always in the basin of at- 
traction of the non-Fermi liquid fixed point and did not 
yet start flowing towards a Fermi liquid one driven by 
the applied field. Since we had to simultaneously cap- 
ture scales that were many orders of magnitude apart 
{h s j vs r) the numerical calculations had to be done us- 
ing "quadruple precision" arithmetics and the algorithms 
had to be designed carefully to assure the precision re- 
quired. Let us also point out that when fitting for in- 
stance xfmp f° r negative values of e (i.e. determining the 
temperature asymptotics of the flavor linear susceptibil- 
ity as the system goes into the overscreened magnetic mo- 
ment regime) , the quantities involved become very small 
and the relative errors grow rapidly. 

Our results show that the infrared physics of the two- 
channel Anderson impurity is described by a line of 
non-Fermi liquid fixed points characterized by a resid- 
ual entropy S® mp = fee In y2 and logarithmic temper- 
ature divergences in the specific heat coefficient (i.e. 



7 = Cj mp /T) and the susceptibilities — all these features 
being already present in the limiting case embodied by 
the two-channel Kondo impurity. What distinguishes the 
different points on the line are, for instance, the chang- 
ing values of the impurity charge valence n c (e — fi) (see 
Fig- EH)- as well as those of the scaling fields and the char- 
acteristic scales. 
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FIG. 3: Schematic renormalization group flow diagram of the 
two-channel Anderson impurity model. The flow connects the 
line of microscopic theories characterized by the ratio e/A 
(in this figure A = T) with the line of infrared fixed points 
parametrized by n c . 



Since the specific heat coefficient and the susceptibili- 
ties have the same singular low-temperature divergences 
along the whole line of fixed points, it is possible to define 
the following two Wilson ratios^ 



Rt'/ = Mm 



S,j j S,J 

Ximp / Xbul] 



T^O Ci m p/Cbu 



(72) 



where s (/) label the magnetic ( quadrupolar) susceptibili- 
ties. This in particular gives for the "flavor (or quadrupo- 
lar) Wilson ratio" : 



i?4 (e - fi) = — lim 



/ 

V. 

A imp 



8/3 



3 t^o C imp /T i 



(A s /A,r 

8/3 
1 + e - 4 ^( £ -^) 



(73) 



with 7 = 0.7838 ± 0.0004. These ratios are another char- 
acteristic that varies along the line of fixed points. In 
particular, we have that i?£ (+oo) = 8/3 in agreement 
with the result for the two-channel quadrupolar Kondo 
model and, as expected, (— oo) = 0, since this is the 
magnetic Kondo limit. By symmetry arguments, we find 
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that Kie-fJ.) = R&(p-f) and R l(0) = 4/3. Fi- 
nally, a simple calculation yields the universal relation: 
Rl (s-fi) + R s w (s-fi) = 8/3. 

The measurement of R^* is a very sensitive way to 
experimentally determine which doublet is the lower one 
in energy, whether the impurity ion ground state is mag- 
netic or quadrupolar. This can be very useful if the sys- 
tem is in the mixed valence region \e — fi\ < T, as was 
thought to be the case in UBei3. This question triggered 
the interest in measuring the nonlinear susceptibility of 
this compound 1 - and its answer is still a source of some 
controversy^ 

E. Nonlinear Magnetic Susceptibility 

The nonlinear magnetic susceptibility X3 measures the 
leading nonlincarity in the magnetization 

M = X B + lx3B 3 + ... (74) 



in the direction of an applied field B. As follows from the 
analysis in Ref. IHM a quadrupolar moment contributes 
nontrivially to X3 via its coupling to the square of the 
magnetic field. In our formalism this coupling adds a 
term aB 2 Jj to the scaling Hamiltonian, where a is a 

constant and J| = (1/2) -ipa^T^ip^^: is the component 
of the quadrupolar moment along the field direction, with 
t z the diagonal SU(2) flavor generator. 

The fact that the quadrupolar moment couples to B 2 
permits the determination of the quadrupolar susceptibil- 
ity from measurements of the (nonlinear) magnetic sus- 
ceptibility. Here we shall explore how the leading con- 
tribution to \3 from the Anderson impurity, X3,imp call 
it, behaves as one goes along the critical line from the 
quadrupolar to the magnetic two-channel Kondo fixed 
points. 

We thus have to consider the extended scaling Hamil- 
tonian 



/oo />oo 
dxJ z s {x) + aB 2 / dxJ z f (x) + subleading terms , (75) 
-oo J — oo 

I 



where we have removed "by hand" terms from the charge 
sector (which do not contribute to X3,imp)- Given 

_ u Jimp I l>7a\ 
X3Srap — ~QB^~ ° 

we now have to carry the linked cluster expansion of 
I 



the free energy to 0(B 4 ), in analogy to the 0(B 2 ) ex- 
pansion in the calculation of the linear susceptibility in 
Sec. (IV. C). To second order in the scaling fields A s and 
Xf we then find that fi mp can receive possible contri- 
butions from 15 different terms generated by H sca i ing in 
Eq. I|75|). These terms have the structure 



a °B A/Ay 3 (n i= oa:---io^( z n j=o,i,... Jo ^/(^+^ , 

where k + £ = 1 or 2, and i + 2j a = 4 (with i = 0,2,4; j = 0,2), and where n m=0 •■• = 1 (m = i,j,k,£). By 
repeated use of the operator product expansions in Eqs. (|r>f|l and (JSSJ, a straightforward analysis reveals that among 
these fifteen terms only one is nonvanishing^ Passing to the finite-T geometry T c , it has the form 



(Jf(Tl,X 1 )Jf(T 2 ,X 2 )GW(T 3 ,0)6W(T i ,0)) T , 

= 9 



1 n 

h(n - 7"3 + ixi)h(T 2 - t 4 + ix 2 ) + h(r 2 - r 3 + ix 2 )h(T! - r 4 + ixi) \ — sin — (t 3 - t 4 ) | _1 , (77) 

J 7T fj 



r 



which is the same connected four-point function that we 
encountered in Eq. Q68f) . but with flavor now taking the 
place of spin. By comparison with the analysis that fol- 
lows in Sec. (IV. C) we immediately conclude that the 
leading contribution to the nonlinear susceptibility has 
the same logarithmic form as that for the linear suscepti- 



bility, but now parametrized by the "quadrupolar" scal- 
ing field Xf and temperature scale Tf. 

X3,i mp = l8X 2 a 2 H^) + 0(X 3 f ). (78) 
Combined with the Bethe Ansatz fit of A/ and Tf from 
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the previous section, Eq. (|78() completely specifies the 
leading scaling term of the nonlinear magnetic suscepti- 
bility. 

V. SUMMARY 

We have carried out in this paper a detailed analy- 
sis of the low-temperature thermodynamics of the two- 
channel Anderson impurity model, identifying and fitting 
a boundary conformal field theory to the Bethe Ansatz 
solution of the model. Combining the methods allowed 
us to completely determine the various parameters and 
scales introduced in the BCFT and, in terms of these, 
give a complete description of the line of fixed points 
characterizing the low-energy physics. 



As we mentioned in the Introduction, having identified 
the scaling operators the asymptotic dynamical proper- 
ties can in principle be calculated. This allows the study 
of resistivities, optical conductivities, and Green's func- 
tions in general. These results will be presented in a 
subsequent work. 
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